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Recently, we have derived a novel and compact expression for how perturbations in the matter 
fields of the cosmological fluid can lead to deviations from the standard Friedmann equations. 
Remarkably, the dissipative damping of velocity perturbations by bulk and shear viscosity in 
the dark sector can modify the expansion history of the universe on arbitrarily large scales. 
In universes in which this effect is sufficiently sizeable, it could account for the acceleration of 
the cosmological expansion. But even if dark matter should be less viscous and if the effect 
would be correspondingly smaller, it may have observable consequences in the era of precision 
cosmology. Here, we review the origin of this backreaction effect and possibilities to constrain 
it further. 


1 Backreaction effects in cosmology 

In the cosmological perturbation theory that is at the basis of the phenomenologically successful 
ACDM model, all fields entering Einstein’s field equations are split into a spatially homoge¬ 
neous and isotropic background part and perturbations around it. For instance, one writes 
9iiu{x) = gfiuir) + 5gfj,u{x) for the metric, e{x) = e(r) + 5e{x) for the energy density, etc., where 
the background fields are defined as spatial averages of the full fields. Up to first order in per¬ 
turbations, the spatially averaged Einstein equations reduce then to the well-known Friedmann 
equations for the background fields. However, Einstein’s field equations are non-linear, and, 
beyond first order in perturbations, a spatial average of Einstein’s equations for the full fields 
differs from Einstein’s field equations for the spatially averaged fields. The perturbations back- 
react, that is, the time evolution of the background fields depends on the time evolution of the 
perturbations. 

The possible sources of such backreaction effects, and the question to what extent they 
may modify Friedmann’s equations have been discussed repeatedly. On the one hand, it has 
been argued that the formation of non-linear structures could affect the average expansion rate 
of the Universe considerably and that this may even account for the accelerated expansion in 
late-time cosmology, for a review see Ref. In sharp contrast to this line of thought, there 
are also arguments that the phenomenological success of the ACDM model is a consequence of 
its mathematical consistency, in the sense that the evolution of the exactly homogeneous and 
isotropic FLRW metric does indeed account with high precision for the evolution of spatially 
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averaged fields in our physically realized Universe with perturbations, for a review see Ref. . 
So far, this discussion focussed mainly on backreactions of metric perturbations, that is, terms 
in Einstein’s field equations that are second order or higher in Instead, we point here 

to backreaction effects that are second order in matter perturbations and that have a physical 
interpretation in terms of dissipative processes. 


2 Backreaction from matter perturbations and cosmological expansion 


Einstein’s equation for an exactly homogeneous and isotropic FLRW universe are the first Fried¬ 
mann equation, that relates e to the Hubble parameter H, and the second Friedmann equation. 
The only additional information contained in the latter is energy conservation for an expanding 
system, 

h+3Hie + p-3CH) = D, (1) 

a 

where D = 0 for an unperturbed FLRW metric, H = a/a^, and the dot denotes a derivative 
w.r.t. conformal time. We recall these elementary facts to highlight two issues that are important 
for the following. First, what enters the evolution of the background helds is not the average 
pressure p, but the effective pressure peff = P — 3(H for which a term proportional to the bulk 
viscosity C times expansion scalar is subtracted. Indeed, bulk viscosity is the only dissipative 
effect that can arise in an exactly homogeneous and isotropic system, and bulk viscous pressure 
is a negative pressure. The logical possibility that bulk viscosity may therefore account for 
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phenomena normally attributed to dark energy has been explored in several studies'^- . Second, 
equation (1) illustrates the well-known fact that Einstein’s field equations include the conditions 
for energy momentum conservation, = 0. In our recent study we started from energy 

conservation. Working in the Landau Frame = eu'^ and keeping shear viscous and bulk 

viscous effects up to first order in a gradient expansion of the fluid dynamic helds, one hnds 
from = 0 

u^dij,e (e -b = 0 , (2) 

where q and C denote the shear and bulk viscosity, respectively, -|- u^u'^ and 
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Forming the spatial average of eq. (2), one can determine the dissipative corrections to the time 
evolution of the average energy density e. In general, metric perturbations enter this spatial 
average via the covariant derivative V^, the projector A^,y and the normalization of the velocity 
held. For instance, in conformal Newtonian gauge with Newtonian potentials T and the 
constraint = —1 implies = {q^qv), where 7 = l/{ay/\ — iP + 2\k -b However, in 

the relevant limit of small gravitational potentials ($ , T <C 1 ), that vary slowly in time, $ ~ ^<I>, 
one hnds that the spatial average of (2) is independent of metric perturbations^. On the other 
hand, there are huid perturbations. In particular, perturbations in the gradients of the velocity 
held enter the terms 0 ^ and in eq. ( 2 ) to second order, and thus they do not vanish 

upon spatial averaging. Since 0^ and are positive semidehnite, one hnds (for small huid 

velocity v) that the spatial average of ( 2 ) is of the form of ( 1 ) with a positive term^ 


D = ^{q diVjdiVj + diVjdjVi - IdiVidjVj ) + ^(C[V • + l{v-V {p- QCH)) . (4) 


Here, the term D gives an explicit backreaction of huid perturbations on the time evolution 
of the average energy density e. In Ref. we have analyzed the effect of this term on the 
cosmological expansion. To this end, we combined eq.(l) with an equation for the scale factor a. 
For the latter, we used the trace of Einstein’s equation since its spatial average does not depend 
on huid velocity perturbations. 


R = S-rrCf^Tlf = —SttGn (e — 3pefr) = . (5) 

Combining both equations, one hnds for the deceleration parameter q = —1 — H/{aH‘^) the 
differential equation^ (we use the simple equation of state pes = we) 
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din a 


+ 2(g - 1) (g - 5(1 -b 3u;)) 
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The fixed point of this equation is accelerating, i.e. q < 0, for > 3 ^ 5 ^ 
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An explicit deviation from Friedmann’s eqnation for the (toy) ansatz D oc H 


The conclusion of the discussion above is that in a universe in which dissipative effects are suffi¬ 
ciently large, D ~ 0(ff^/GN), dissipation induces accelerated cosmological expansion. Clearly, 
in a universe in which dissipative effects are smaller, other physical effects need to be invoked 
to account for an observed accelerated expansion. However, even for small but non-vanishing 
backreaction D, the cosmic expansion scale factor a will not obey the exact Friedmann equation. 
We emphasize this point since it is a key assumption of the ACDM cosmological concordance 
modeP that a does satisfy the Friedmann equation; this assumption underlies for instance the 
definition of cosmological distance measures and it is poorly tested at best. Given the precision 
of modern cosmology, it is thus interesting to explore the sensitivity of distance measures on D 
by analyzing deviations from Friedmann’s equation induced by D. 

The term D is calculable once the matter content of the Universe and the dissipative proper¬ 
ties of dark matter are specified. D affects the time evolution of the background fields in eq. (1), 
while the cosmological perturbations needed to determine the spatial average D according to 
eq. (4) are evolved on top of these background helds. This complicates a dynamical determi¬ 
nation of D. In general, since perturbations evolve, D will have a non-trivial dependence on 
redshift, D = D{z), where a = We plan to calculate this evolution of D in future work for 

different viscous properties of the dark sector ‘. For the purpose of these proceedings, we just 
examine heuristically one type of z-dependence of D: 

Let us consider for a moment a hypothetical universe composed of radiation and viscous 
dark matter, but without dark energy. What would we have to choose for the dissipative 
backreaction term D of that hypothetical universe, if we required that its deceleration pa¬ 
rameter matches that of the phenomenologically successful ACDM universe? To answer this 
question, we write the scale dependence of the Hubble parameter of the ACDM universe, 
H{z) = Hoa/Ha -|- we derive from it the explicit form for the de¬ 

celeration parameter q{z) in the ACDM universe, we insert it into the left hand side of equation 
( 6 ) and we read off an expression for D. For the case of a pressure-less fluid = 0, we find 
D = Therefore, if one would aim at exploring the possibility that dissipative 

properties of dark matter can account for phenomena normally attributed to dark energy, one 
would seek dark matter with viscous properties that lead to a backreaction D(z) ~ const ■ H[z). 

We caution that by making the ad hoc ansatz D{z) = const- H{z), we did not check whether 
this particular scale dependence can arise from a dynamical evolution of perturbations in a fluid 
with ’realistic’ viscous properties. In general, both the size and the scale dependence of D{z) 
depend on the viscous properties of dark matter. There is no a priori reason that the scale 
dependence satisfies the simple ansatz D{z) = const • H{z). For a more general form of D{z), 
the solution of the evolution equations for the background field e will show deviations from 
Friedmann’s equation that involve an integration over the evolution of D. However, for the 
ansatz D{z) = const • H[z), a linear superposition of eqs. (1) and (5) results in 
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Here, the relation between the Hubble parameter and the total energy density e and total 
effective pressure differ from the standard Friedmann equations by the terms proportional to D. 
There is a one-to-one relation between this type of evolution and a ACDM model with energy 
density e = e — ^ and dark energy density y|^. More precisely, in terms of e, the evolution 
equation (1) for energy density, and the trace (5) of Einstein’s equations can be written as 


1 • SttGn D 

-e + 3H (e + Pes) = 0 , and R-\ - 
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3H 


SttGn (e - 3pes) . 


( 8 ) 










This is of the same form as the corresponding equations in the ACDM model if one specializes 
to a pressureless fluid peff = 0 and if one identifies e with the energy density in the ACDM model 
and the cosmological constant with 


3H 


(9) 


4 Outlook 

Rather than repeating or anticipating work published in refereed journals, we have made use 
of this proceedings article to share some simple but heuristic considerations of how dissipative 
backreaction of matter perturbations may lead to deviations from Friedmann’s equation for 
the scale parameter. Clearly, a more rigorous analysis of these deviations is possible: one 
can calculate the size and scale dependence of D[z) as a function of the viscous properties of 
dark matter. This calculation could contribute to a critical appraisal of cosmological distance 
measures. This line of thought will be presented elsewhere. 

We finally remark on a crucial point that was discussed more prominently in the oral presen¬ 
tation of our work in Moriond. The evolution of perturbations and the growth of cosmological 
structure depends on dissipative properties, and information about structure formation is known 
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to constrain possible dissipative phenomena"^’ . To first order in cosmological perturbations, 
one finds from Einstein equations that the evolution equation of the density contrast 5p = 5e/e 
does not depend explicitly on viscous contributions, while the evolution of fluid velocity per¬ 
turbations is attenuated by a scale-dependent viscous term. As a consequence, one expects 
that unlike finite pressure that counteracts the growth of density contrast, finite viscosity will 
not counteract this growth but it will slow down the growth of velocity gradients and it may 
thus delay structure formation. We are currently evaluating ‘ how the constraints on viscous 
properties arising from structure formation translate into constraints on the backreaction D. 
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